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Abstract 

This  paper  describes  a  newly  devised  autosomal  short  tandem  repeat  (STR)  multiplex  polymerase  chain  reaction  (PCR)  systems  for  10 
loci  (D1S1656,  D2S1353,  D8S1132,  D12S1090,  D14S608,  D18S535,  D19S253,  D20S480,  D21S226,  and  D22S689)  unlinked  to  the  core 
STR  loci  (non-CODIS  loci).  Of  252  samples  taken  from  the  Japanese  population,  PCR  products  ranged  in  length  from  107  bp  to  319  bp. 
No  significant  deviations  from  Hardy-Weinberg  equilibrium  were  observed  at  any  of  the  10  loci.  The  accumulated  power  of  discrimi¬ 
nation  and  power  of  exclusion  for  the  10  loci  were  0.999999999998  and  0.99991,  respectively.  We  conclude  that  the  present  multiplex 
system  for  the  10  non-CODIS  loci  represents  a  powerful  tool  for  forensic  applications. 

©  2008  Elsevier  Ltd  and  FFLM.  All  rights  reserved. 
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Introduction 

Analyses  of  short  tandem  repeat  (STR)  markers  are  an 
essential  component  of  forensic  genetics.  Since  the  13  core 
combined  DNA  index  system  (CODIS)  STR  loci  were  used 
for  a  DNA  database  of  felons,  several  commercially  avail¬ 
able  multiplex  polymerase  chain  reaction  (PCR)  kits  have 
been  developed  for  analysis  of  core  STR  loci  and  have 
now  entered  wide  use.  On  the  other  hand,  many  highly 
polymorphic  STR  loci  exist,  which  are  unlinked  to  the  core 
STR  loci  (non-CODIS  loci).  Analysis  of  such  polymorphic 
non-CODIS  loci  in  addition  to  the  core  CODIS  loci  may 
increase  the  power  of  exclusion  in  complex  kinship  testing 
involving  deficient  paternity  cases.  The  present  study 
sought  to  develop  a  new  multiplex  PCR  system  for  simul¬ 
taneous  typing  of  highly  polymorphic  non-CODIS  loci. 
We  selected  10  polymorphic  non-CODIS  loci  (D1S1656, 
D2S1353,  D8S1132,  D12S1090,  D14S608,  D18S535, 
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D19S253,  D20S480,  D21S226,  and  D22S689)  from  the  lit¬ 
erature,  and  surveyed  their  efficacy  for  forensic  practice. 

Materials  and  methods 

Informed  consent  was  obtained  and  blood  samples  col¬ 
lected  from  252  healthy  unrelated  individuals  residing  in 
Japan,  and  additionally  24  daughters  and  21  sons  (45  fam¬ 
ily  trios).  Parents  of  the  children  were  included  among  the 
252  subjects,  and  the  paternal  and  maternal  relationships 
of  the  families  were  confirmed  by  autosomal  STRs  geno- 
typing  using  the  AmpflSTR  Identifier  Kit  (Applied  Biosys¬ 
tems,  Foster  City,  CA).  DNA  was  extracted  from  whole 
blood  using  the  nucleic  acid  isolation  system  Quick 
Gene-800  (FUJIFILM,  Tokyo,  Japan)  according  to  the 
protocol  described. 

Multiplex  PCR  for  amplification  of  the  10  non-CODIS 
STR  loci  was  performed  using  the  fluorescent  dye  labeled 
primer  sets  listed  in  Table  1.  PCR  reactions  were  performed 
in  a  total  volume  of  15  pi  containing  of  2.0  ng  of  genomic 
DNA,  lx  GeneAmp  PCR  buffer,  1.5  mM  of  MgCl2, 
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Table  1 

Primer  sequences  used  in  this  study 
STR  locus  Primer  sequences 

Forward  primer 

D 1 S 1 656  5'-GCTCAAGGGTCAACTGTAT-3' 

D2S 1353  5'-AGCCAGGGACATTGCTTAA-3' 

D8S 1 1 32  5'  -GGCT  AGG  AA  AGGTT  AGT  GGC-3' 

D 1 2S 1 090  5'-AAGCCCAAAGATGTAAGGCT-3' 

D14S608  5'  -T  AAAGGTTT  AT  CCATGCTGT  AGC-3' 

D18S535  5'-TCATGTGACAAAAGCCACAC-3' 

D19S253  5'  -  AT  AG  AC  AGAC  AG  ACGG  ACT  G-  3' 

D20S480  5'-GTGGTGAACACCAATAAATGG-3' 

D21S226  5'-GGAAACCACTCTAAGACATA-3' 

D22S689  5'-TATGTACAGACCTGCAACTTGC-3' 


Reverse  primer11 

5'-AAATAGAATCACTAGGGAA-3' 

5'-CCATTAGAGCAGGATTTGT-3' 

5'-CCCTCTCTCTTTCGAGCAAT-3' 

5'-ACCAACCTAGGAAACACAGT-3' 

5'- ACGT  GGT  AC  AGGT  AGAT  A  AAT  GG-3' 

5'- AG  AC  AGAAAT  AT  AG  AT  G  AG  AAT  GC  A-3' 
5'-GGG  AGT  GGAG  ATT  ACCCCT  -  3' 
5'-AAGCAAATAAAACCAATAAACTCG-3' 

5'-CCTGCCTGCCTATCTATCTG-3' 


Dye  label 

Reference 

VIC 

Newly  designed’ 

PET 

Newly  designed’ 

NED 

GDBC 

VIC 

7 

PET 

GDB 

FAM 

GDB 

FAM 

GDB 

NED 

GDB 

PET 

Newly  designed1 

NED 

GDB 

a  A  particular  sequence  of  seven  bases,  a  tailed  primer  (Applied  Biosystems),  was  added  to  the  5'  end  of  all  reverse  primers  to  promote  adenylation. 
b  Newly  designed  primer  set  in  this  study. 

0  Genome  Database  (http://www.gdb.org). 

d  Newly  designed  reverse  primer  in  this  study.  Forward  primer  is  shown  in  the  GDB. 


200  pM  of  each  deoxyribonucleotide  triphosphate  (dNTP, 
GeneAmp  dNTP  MIX),  each  primer  set,  and  1.25  U  of 
AmpliTaq  Gold  DNA  polymerase  (Applied  Biosystems). 
The  multiplex  PCR  mixture  contained  0.25  pM  primer  sets 
of  D18S535  (fluorescent  dye  labeled  with  6FAM),  0.12  pM 
D19S253  (6FAM),  0.5  pM  D1S1656  (VIC),  0.25  pM  D12- 
S1090  (VIC),  0.25  pM  D8S1132  (NED),  0.65  pM  D22- 
S689  (NED),  1.0  pM  D20S480  (NED),  0.95  pM  D21S226 
(PET),  1.0  pM  D2S1353  (PET),  and  0.4  pM  D14S608 
(PET).  The  cycling  programs  consisted  of  pre-denaturation 
at  95  °C  for  1 1  min,  followed  by  30  cycles  of  denaturing  at 
94  °C  for  30  s,  annealing  at  56  °C  for  30  s,  extension  at 
72  °C  for  45  s,  and  a  final  extension  at  70  °C  for  60  min  with 
a  GeneAmp  9700  (Applied  Biosystems)  in  9600  mode.  Elec¬ 
trophoresis  was  performed  with  an  ABI  3100-Avant 
Genetic  Analyzer  (Applied  Biosystems)  in  which  1 .0  pi  of 
PCR  product  was  mixed  with  15  pi  Hi-Di  formamide  and 
0.2  pi  of  the  GeneScan-500LIZ  size  standard.  GeneMapper 
ID  v3.2  Software  was  used  to  analyze  the  data.  Allelic  typ¬ 
ing  was  based  on  sequenced  allelic  ladders.  All  fragments 
contained  in  allelic  ladder  were  sequenced  to  calibrate  the 
number  of  repeats.  Non-labeled  primer  sets,  which  had 
the  same  sequences  as  given  in  Table  1,  were  used  for 


sequencing.  The  target  alleles,  re-amplified  after  extraction 
using  the  crush  and  soak  method,1  were  used  as  templates. 
Sequencing  reactions  were  performed  with  the  same  primers 
using  the  BigDye  vl .  1  Ready  Reaction  Kit  (Applied  Biosys¬ 
tems).  Following  electrophoresis  with  the  ABI  3100-Avant 
Genetic  Analyzer  (Applied  Biosystems),  sequencing  was 
determined  with  sequencing  analysis  software  5.1.1  (Applied 
Biosystems).  Genotypes  for  252  Japanese  individuals  and 
commercially  available  cell  lines,  including  9947A  (Applied 
Biosystems),  K562  and  9948  DNA  (Promega,  Madison, 
WI)  were  typed  using  the  allelic  ladder.  To  determine  the 
minimum  quantity  of  DNA  required  to  achieve  reliable 
results  [peak  highs  at  all  loci  exceeded  100  relative  fluores¬ 
cence  units  (RFU)]  with  the  multiplex  systems,  we  used  ali¬ 
quots  of  the  cell  line  9947A  DNA  (50,  100,  200,  300,  400, 
500,  750  pg,  1,  5,  and  10  ng). 

The  allele  frequency  at  each  locus  was  calculated  from 
the  numbers  of  alleles  in  the  252  Japanese  individuals.  Var¬ 
ious  forensic  parameters  were  calculated  using  the  Power- 
Stats  vl.2  software  package  (Promega).  Hardy-Weinberg 
equilibrium  and  population  differentiation  were  deter¬ 
mined  by  an  exact  test  using  the  GENEPOP  (version  3.4) 
software  package  designed  by  Raymond  and  Rousset.2 


Table  2 

Information  on  the  STR  loci  in  this  study  and  a  result  of  genotype  of  cell  line  DNA 


STR  locus 

Product  length  (bp)a 

Allele  range 

Number  of  alleles 

Genotype  of  cell  line  DNA 

K562  9947A 

9948 

Chromosomal  location 

D1S1656 

129-164 

11-19.3 

13 

15,  16 

18.3,  18.3 

14,  17 

lq42 

D2S1353 

150-180 

8-18 

11 

10,  12 

12,  14 

12,  14 

2q24 

D8S1132 

142-178 

16-25 

10 

20,  24 

19,  21 

20,  24 

8q23 

D12S1090 

213-297 

8-29 

20 

21,  25 

19,  20 

21.  22 

12q  12 

D14S608 

201-233 

6-14 

9 

7,  7 

7.  11 

12,  12 

14q  12 

D18S535 

138-174 

9-18 

10 

13,  14 

13,  15 

13,  15 

1 8q  12 

D19S253 

213-249 

6-15 

10 

11,  12 

7,  8 

11,  12 

19pl  3. 1 

D20S480 

279-319 

10-20 

11 

16,  17 

15,  16 

16,  17 

20ql3.2 

D21S226 

107-123 

7-11 

5 

10,  10 

10,  10 

9.  10 

2 1  q2 1 

D22S689 

202-234 

12-20 

9 

16,  16 

16,  18 

17,  17 

22ql2 

a  Product  length  includes  seven  bases  sequence  (tailed  primer)  and  an  adenine  due  to  adenylation. 
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Results  and  discussion 

We  reviewed  the  literature  and  selected  10  polymorphic 
non-CODlS  STR  loci  (D1S1656,3  D2S1353,4  D8S1132,5 
D12S1090,6’7  D14S608,8  D18S535,9  D19S253,10  D20- 
S480,11'12  D21S226,12  and  D22S68912)  located  on  discrete 
chromosomes  (Table  2).  D1S1656,  D14S608,  D19S253, 
D20S480,  and  D22S689  are  found  on  different  chromo¬ 


somes  from  the  CODIS  core  loci.  Although  D2S1353  and 
D12S1090  are  found  on  chromosomes  forming  part  of 
the  CODIS  core  loci  (TPOX  and  vWA,  respectively),  the 
loci  are  located  on  different  arms,  and  they  are  separated 
by  a  considerable  distance.  D8S1132,  D18S535,  and 
D21S226  are  found  on  the  same  arms  as  the  CODIS  loci 
(D8S1179,  D18S51,  and  D21S11,  respectively),  but  these 
loci  are  more  than  10  Mb  apart  from  one  another.  Genetic 


9947 A  j  lOmultiplex  |  j  g 

1  D21S226  1  I  D2S1353  I  I  D14S608  I 

100  140  180  220  260  300 


14 


Fig.  1.  Electropherogram  of  DNA  typing  result  of  control  DNA  9947A  using  the  present  multiplex  system  for  assay  with  10  non-CODIS  loci. 
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linkage  disequilibrium  (LD)  in  the  human  genome  is  gener¬ 
ally  reported  at  distances  of  50-100  kb. 13-15  Thus,  LD 
between  the  current  non-CODIS  eight  loci  and  the  CODIS 
core  loci  may  be  excluded  from  consideration.  Analysis  of 
extended  STR  loci  using  both  the  present  multiplex  system 
and  commercial  multiplex  kits  that  include  the  CODIS  core 
loci  have  been  demonstrated  to  be  remarkably  effective  in 
cases  involving  difficult  kinship  testing.  The  10  non-CODIS 
loci  examined  in  this  study  were  analyzed  using  a  newly 
developed  multiplex  PCR  system  that  included  10  primer 
sets  customized  by  Applied  Biosystems.  A  7-base  sequence 
(tailed  primer,  Applied  Biosystems)  was  appended  to  the  5' 
end  of  the  all  reverse  primers  to  promote  adenylation. 16-18 
PCR  amplification  and  allele  typing  with  the  system  was 
successful  for  all  10  loci.  Fig.  1  shows  an  example  of  the 
typing  of  a  commercially  available  cell  line  9947A  DNA. 
Very  minor  artifact  peaks  were  detected,  which  had  no 
effect  on  allele  typing.  Following  analysis  of  the  sensitivity 


range  of  the  DNA  detection  limit,  allele  typing  at  all  10  loci 
proved  successful  when  using  at  least  400  pg  of  template 
DNA  with  peak  heights  exceeded  100  RFU.  In  an  analysis 
with  200  and  300  pg  of  template  DNA,  typing  at  all  10  loci 
was  successful  due  to  an  increase  in  the  amount  of  Taq 
polymerase  (2.5  U).  The  amplification  products  for  the  10 
loci  from  252  Japanese  individuals  ranged  in  length  from 
107  bp  for  D21S226  to  319  bp  for  D20S480.  Numbers  of 
alleles  ranged  from  5  alleles  for  D21S226  to  20  alleles  for 
D12S1090.  Table  2  provides  information  on  the  10  loci 
and  results  for  the  genotype  of  a  cell  line  DNA,  while  the 
whole  genotyping  data  of  252  Japanese  individuals  are 
shown  in  Supplementary  Table  1 .  Previous  papers  have  used 
two  styles  of  allelic  nomenclature  for  D20S480.  One  paper 
refers  to  the  repeat  motif  of  allele  (n)  as  (TCTA)n,u 
while  the  other  refers  to  (TATC)n_2ATC(TATC)2.12  The 
present  study  followed  the  latter  convention  for  nomen¬ 
clature. 


Table  3 

Allele  frequencies  and  statistical  parameters  of  the  10  loci  in  a  Japanese  population  (JV  =  252) 


Allele 

D1S1656 

D2S1353 

D8S1132 

D12S1090 

D14S608 

D18S535 

D19S253 

D20S480 

D21S226 

D22S689 

6 

0.040 

0.002 

7 

0.192 

0.153 

0.008 

8 

0.002 

0.004 

0.026 

0.046 

0.087 

9 

0.002 

0.008 

0.139 

0.234 

0.010 

0.143 

10 

0.050 

0.052 

0.216 

0.022 

0.028 

0.002 

0.746 

11 

0.028 

0.018 

0.133 

0.206 

0.008 

0.121 

0.002 

0.016 

12 

0.052 

0.143 

0.117 

0.135 

0.131 

0.333 

0.079 

0.026 

13 

0.109 

0.290 

0.038 

0.044 

0.270 

0.246 

0.052 

0.006 

14 

0.067 

0.139 

0.020 

0.002 

0.238 

0.058 

0.117 

0.012 

15 

0.234 

0.175 

0.077 

0.004 

0.222 

0.133 

16 

0.284 

0.083 

0.004 

0.016 

0.218 

0.240 

16.3 

0.012 

17 

0.087 

0.095 

0.034 

0.004 

0.002 

0.228 

0.363 

17.3 

0.063 

18 

0.014 

0.004 

0.212 

0.004 

0.002 

0.069 

0.161 

18.3 

0.040 

19 

0.004 

0.196 

0.008 

0.008 

0.052 

19.3 

0.006 

20 

0.165 

0.002 

0.002 

0.008 

21 

0.117 

0.006 

22 

0.143 

0.014 

23 

0.099 

0.032 

24 

0.028 

0.069 

25 

0.002 

0.147 

26 

0.157 

27 

0.141 

28 

0.038 

29 

0.008 

OH11 

0.825 

0.821 

0.845 

0.865 

0.877 

0.798 

0.766 

0.766 

0.413 

0.778 

PDb 

0.952 

0.946 

0.952 

0.976 

0.946 

0.922 

0.916 

0.945 

0.623 

0.908 

PEC 

0.647 

0.639 

0.686 

0.725 

0.749 

0.595 

0.537 

0.537 

0.122 

0.558 

PICd 

0.81 

0.81 

0.82 

0.88 

0.81 

0.76 

0.75 

0.80 

0.38 

0.73 

TPIe 

2.86 

2.80 

3.23 

3.71 

4.06 

2.47 

2.14 

2.14 

0.85 

2.25 

P1 

0.5799 

0.1414 

0.3366 

0.5896 

0.8610 

0.4140 

0.0660 

0.1540 

0.5509 

0.8754 

a  Observed  heterozygosity. 
b  Power  of  discrimination. 

0  Power  of  exclusion. 
d  Polymorphism  information  content. 
e  Typical  paternity  index. 

1  P-values  of  the  exact  tests  for  Hardy-Weinberg  equilibrium. 
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Table  3  shows  the  observed  allele  frequencies  and  foren¬ 
sic  statistical  parameters  for  the  10  loci  in  the  Japanese 
population.  No  significant  deviations  from  Hardy-Wein- 
berg  equilibrium  were  observed  in  any  of  the  10  loci 
(P  >  0.05).  Heterozygosity  values  exceeded  0.76,  except 
for  D21S226  (0.413).  The  power  of  discrimination  values 
ranged  from  0.976  for  D12S1090  to  0.908  for  D22S689, 
except  for  D21S226  (0.623).  All  forensic  parameters  for 
D21S226  were  lowest  among  the  10  loci,  while  the  other 

9  loci  exhibited  high  polymorphism.  The  accumulated 
power  of  discrimination  and  power  of  exclusion  for  the 

10  loci  were  0.999999999998  and  0.99991,  respectively. 
No  mutations  were  detected  in  the  kinship  cases  investi¬ 
gated  (45  family  trios).  Given  the  infrequent  meioses 
observed  in  the  study,  estimating  a  reliable  mutation  rate 
would  require  further  investigation. 

The  previous  literature  report  several  earlier  population 
studies  of  the  present  10  non-CODIS  loci  from  Asian  pop¬ 
ulations.  Single  locus  comparisons  from  this  published 
data  indicate  that  the  distribution  of  D19S253  alleles  is 
similar  to  those  in  Korean  and  Thai  populations19'20  ( P 
values  =  0.65  and  0.30,  respectively).  In  D14S608,  the  Kor¬ 
ean  population8  differed  significantly  from  the  present  Jap¬ 
anese  population  (P  <  0.00001).  Comparisons  of  available 
Chinese  published  data21  22  indicate  that  distributions  of 
D2S1353  and  D18S535  alleles  also  differ  significantly  from 
the  Japanese  data  ( P  <  0.00001).  On  the  other  hand,  pub¬ 
lished  data  on  the  German  population  exist  for  the  eight 
loci,6'12’23’24  except  for  D2S1353  and  D19S253  and  the  dis¬ 
tribution  of  alleles  for  D8S1132  alone  showed  no  differ¬ 
ences  between  German  and  Japanese  populations 
(P  =  0.32).  However,  significant  differences  between  Ger¬ 
man  and  Japanese  populations  were  found  at  the  other  loci 
(P<  0.00001). 

We  conclude  that  analysis  of  extended  STR  loci  using 
the  present  multiplex  PCR  system  for  ten  non-CODIS  loci, 
when  added  to  analysis  of  core  CODIS  loci,  can  be  highly 
effective  for  difficult  kinship  testing,  including  deficient 
paternity  cases. 

Appendix  A.  Supplementary  material 

Supplementary  material  associated  with  this  article  can 
be  found,  in  the  online  version,  at  doi:10.1016/j.jflm.2008. 
04.001. 
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